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The four standard quantum optics models such as Rabi, Dicke, Jaynes-Cummings ( JC ) and 
Tavis-Cummings (TC) model were proposed by the old generation of great physicists many decades 
ago. Despite their relative simple forms and many previous theoretical works, their solutions at 
a finite N, especially inside the superradiant regime, remain unknown. In this work, we address 
this outstanding problem by using the 1/J expansion and exact diagonization to study the [7(1)/Z 2 
Dicke model at a finite N. This model includes the four standard quantum optics model as its 
various special limits. The 1/J expansions is complementary to the strong coupling expansion used 
by the authors in arXiv: 1512.08581 to study the same model in its dual Z2/U{1) representation. 

We identify 3 regimes of the system’s energy levels: the normal, U{1) and quantum tunneling (QT) 
regime. The system’s energy levels are grouped into doublets which consist of scattering states and 
Schrodinger Cats with even ( e ) and odd ( o ) parities in the U{1) and quantum tunneling (QT) 
regime respectively. In the QT regime, by the WKB method, we find the emergencies of bound states 
one by one as the interaction strength increases, then investigate a new class of quantum tunneling 
processes through the instantons between the two bound states in the compact photon phase. It is 
the Berry phase interference effects in the instanton tunneling event which leads to Schrodinger Cats 
oscillating with even and odd parities in both ground and higher energy bound states. We map out 
the energy level evolution from the U{1) to the QT regime and also discuss some duality relations 
between the energy levels in the two regimes. We also compute the photon correlation functions, 
squeezing spectrum, number correlation functions in both regimes which can be measured by various 
experimental techniques. The combinations of the results achieved here by 1/J expansion and 
those in arXiv: 1512.08581 by strong coupling method lead to rather complete understandings of the 
U{l)/Z 2 Dicke model at a finite N and any anisotropy parameter /3. Experimental realizations and 
detections are presented. Connections with past works and future perspectives are also discussed. 


I. INTRODUCTION 

Quantum optics is a subject to describe the atom- 
photon interactionsiii^. The history of quantum optics 
can be best followed by looking at the evolution of quan¬ 
tum optics models to study such interactions. In the Rabi 
model^, a single mode photon interacts with a two level 
atom with equal rotating wave (RW) and counter rotat¬ 
ing wave (CRW) strength. To study possible many body 
effects such as ” optical bombs”, a single two level atoms 
in the Rabi model was extended to an assembly of N 
two level atoms in the Dicke model^. When the coupling 
strength is well below the transition frequency, the CRW 
term in the Rabi model is effectively much smaller than 
that of RW, so it was dropped in the Jaynes-Cummings ( 
JC ) model^. Similar to the generalization from the Rabi 
to the Dicke model, the single two level atom in the JC 
model was extended to an assembly of N two level atoms 
in the Tavis-Cummings (TC) model^. 

The importance of the 4 standard quantum optics 
model at a finite N in quantum and non-linear optics 
ranks the same as the bosonic or fermionic Hubbard mod¬ 
els and Heisenberg models in strongly correlated electron 
systems and the Ising models in Statistical mechanicsii^. 
There have been extensive theoretical investigations on 


the solutions of the four standard quantum optics model. 
Most of the theoretical works focused on the thermody¬ 
namic limit N ^ 00. The TC model was studied at 
N ^ 00 in^rJ^. A normal to a superradiant phase tran¬ 
sition was found and the zero mode due to the broken 
U{ 1 ) symmetry identified in the superradiant phase. The 
Dicke model at iV —>■ 00 was investigated ini^. A su¬ 
perradiant phase transition with the broken Z2 symme¬ 
try was found and two gapped modes due to the broken 
Z2 symmetry identified in the superradiant phase. How¬ 
ever, there were only very limited works at a finite N. 
The Exact Diagonization (ED) ini^ shows that the level 
statistics in a give parity sector changes from the Pois- 
sonian distribution in the normal phase to the Wigner- 
Dyson in the superradiant phase at any finite N. For 
the Dicke model, the ground state photon number at the 
normal to the superradiant quantum critical point QCP 
was foundi^d^ to scale as {riph) ^ which is a direct 

consequence of finite size scaling near a QCP with infi¬ 
nite coordination numbersidiil. There are also formally 
’’exact” Bethe Ansatz-like solution for the integrable TC 
model at a finite N—. Recently, a formal ’’exact” solu¬ 
tion was found even for the non-integrable Dicke model 
at A^ = 1 ( Rabi model Unfortunately, these ’’exact” 
solutions are essentially useless in extracting any physical 
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FIG. 1. (a) Cold or thermal atoms in a trap embedded 

in a cavity in the strong coupling regime with a transverse 
pumping2i^. The probe detects the Fluorescence spectrum 
of the cavity leaking photons, (b) = 2 ~ 9 superconducting 

qubits are placed on the one or more anti-nodes of a circuit 
QED resonatos^^^. 


phenomen a^^i^^ . 

It is convenient to classify the four well known quan¬ 
tum optics models from a simple symmetry point of view: 
the TC and Dicke model as U{1) and Z 2 Dicke model re¬ 
spectively, while JC and Rabi model are just as the iV = 1 
version of the two^“— . In fact, as stressed in^i, there are 
also two different representations on all of the 4 models: 
N— representation with N independent two level atoms 
with a large Hilbert space 2^ and spin J = N/2 represen¬ 
tation with a smaller Hilbert space TV -|- 1. The relations 
between the two representations were clarified in^i. The 
dramatic finite size effects such as the Berry phase effects, 
Goldstone and Higgs modes on t/(l) Dicke model were 
thoroughly discussed in both N— representation in22 by 
1/A^ expansion and spin-J = N/2 representational by 
1/J representation. Remarkably, we find nearly perfect 
agreements between the results achieved by 1 /N and 1/J 
and the ED studies even when N gets to its lowest value 
TV = 1. The effects of a small CRW term near the C/(l) 
Dicke model limit was also studied inll. In view of the 
tremendous success of the 1 /TV and 1/J expansion in 
studying many strongly correlated electron systems!!^, 
particularly in the t/(l) Dicke model achieved inlS^— , it 
is natural to apply them to study the Z 2 Dicke model as 
originally planned inlS. 

Due to recent tremendous advances in technologies, 
the 4 standard quantum optics models were successfully 
achieved in at least two experimental systems (1) with a 
BEG of TV ^ 10® atoms inside an ultrahigh-finesse 
optical cavity2^“— and (2) superconducting qubits inside 
a microwave circuit cavity^^— or quantum dots inside 
a semi-conductor microcavity^. The superradiant phase 
in the Dicke model was also realized in system (1) with 
the help of transverse pumping ( FigH^ ^ould 

also be realized ’’spontaneously ” in the system (2) with¬ 
out external pumping2i. Indeed, by enhancing the in¬ 
ductive coupling of a flux qubit to a transmission line 
resonator, a remarkable ultra-strong coupling with indi¬ 
vidual g ^ 0.12uJa was realized in a circuit QED systent^^. 
However, in such a ultra-strong coupling regime, the 
system is described well neither by the TC model nor 
the Dicke model, but a combination of the two EqnlT] 
with unequal RW and CRW strength dubbed as U{l)/Z 2 
Dicke model in2i. It was also proposed in^ that in the 


thermal or cold atom experiments, the strengths of g 
and g' can be tuned separately by using circularly po¬ 
larized pump beams in a ring cavity. Indeed, based on 
the scheme, a recent experiments^ realized the U{l)/Z 2 
Dicke model with continuously tunable g and g'. As ar¬ 
gued inSSiSij with only a few TV = 2 ~ 9 qubits em¬ 
bedded in system (2), the hnite size effects may become 
important and experimentally observable. With the re¬ 
cent advances of manipulating a few to a few hundreds 
of cold atoms inside an optical cavity in system 
the finite size effects may also become important and ex¬ 
perimentally observable in near future experiments on 
system (1). As advocated inSi, the Hamiltonian Eqn[T] 
with independent g and g' is the generic Hamiltonian 
describing various experimental systems under the two 
atomic levels and a single photon mode approximation. 
In^i, by the 1/J expansion, we focused on the U[\)/Z 2 
Dicke model EqnlT] near the U{1) limit ( namely, with a 
small anisotropy parameter g'/g = /3 <C 1 and not too 
far from the critical strength gc ) at any finite TV. In 
a very recent preprint^S^ by the strong coupling expan¬ 
sion and the ED, the authors studied the 17(l)/^2 Dicke 
model in its dual presentation starting from from the Z 2 
limit /3 = 1. Here, by the 1/J expansion and ED, we 
will study the U{l)/Z 2 Dicke model Eqn[T] starting from 
from the H(l) limit j3 — Q which is complementary to 
the strong coupling expansion iii^. The combinations 
of both approaches will lead to rather complete under¬ 
standings of the U{l)/Z 2 Dicke model EqnlT] in the full 
range of 0 < /3 < 1. 

In this paper, we study novel quantum phenomena in 
the U{\)/Z 2 Dicke model EqnlT] in its spin J = TV /2 rep¬ 
resentation at a finite TV, any interaction strength g and 
anisotropy parameter 0 < ,5 < 1 by the 1/J expansion^! 
and the ED ^® As a fixed jS, as the g increases, we iden¬ 
tify 3 crossover regimes: the normal, 17(1) and the quan¬ 
tum tunneling (QT) regime Figl^ The super-radiant 
regime at TV = 00 splits into the two regimes at a finite 
TV: the 17(1) and quantum tunneling (QT) regime. In 
the 17(1) regime, we perform a (non-)degenerate pertur¬ 
bation to evaluate the energy spectrum. It is the Berry 
phase which leads to the level crossings between the even 
and odd parity, therefore the alternating parities on the 
ground state and excited states. In the QT regime, by the 
WKB method, we find the emergencies of bound states 
one by one as the interaction strength increases, then 
investigate a new class of quantum tunneling processes 
through the instantons between the two bound states in 
the compact photon phase. It is the Berry phase in¬ 
terference effects in the instanton tunneling event which 
leads to Schrodinger Cats oscillating with even and odd 
parities in both ground and higher energy bound states. 
We map out the energy level evolution from the 17(1) to 
the QT regime. In the 17(1) regime, the doublets con¬ 
sist of scattering states organized as (e, o),(o, e). ( or 

(o, e),(e,o). ) pattern. While in the QT regime, the 

doublets consist of bound states ( or Schrodinger Cats ) 
organized as (e, o), (e, o). ( or (o, e), (o, e). ) pattern. 
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FIG. 2. The three regimes of the U{1)IZ2 Dicke model at 
a finite N as the coupling g increases at a fixed anisotropy 
parameter 0 < g'/g = j3 < 1. The normal regime, the 17 ( 1 ) 
regime and the quantum tunneling (QT) regime. The quan¬ 
tum critical point ( QCP ) at = oo is at gc = (l+P). 

The superradiant regime at = oo splits into the 17 ( 1 ) regime 
and the QT regime at a finite N. The dynamic variables in 
the three regimes are: a,b in the normal regime, (5pa,e®®“ in 
the 17 ( 1 ) regime and Spa, Oa,Tz in the QT regime. As shown in 
the text, there exists some duality relations between the 17 ( 1 ) 
regime and the QT regime. The strong coupling expansion 
ir»2^ works very well in the QT regime and also 17 ( 1 ) regime 
not too close to the QCP. The 17 ( 1 ) regime may disappear if 
P gets too close to 1. 

There are also some sort of duality relations between the 
two regimes in both Hamiltonian and spectrum. We com¬ 
pute the photon correlation functions, squeezing spec¬ 
trum and number correlation functions in both regimes 
which can be detected by Fluorescence spectrum, phase 
sensitive homodyne detection and Hanbury-Brown-Twiss 
(HBT) type of experiments respectively^. When com¬ 
paring the results with those achieved from the strong 
coupling expansion in^, we find nearly perfect agree¬ 
ments among the 1/J expansion, the strong coupling 
expansion and the ED not only in the QT regime, but 
also in the 17(1) regime not too close to the QCP at 
N = oo. The combination of the three methods lead to 
rather complete physical pictures in the whole crossover 
regimes from the 17(1) Dicke to the Z 2 Dicke model in 
Fig|2] Experimental realizations, especially the prepa¬ 
rations and detections of the Schrodinger Cats in both 
experimental systems are discussed. Connections with 
the previous works are speculated and future perspec¬ 
tives are outlined. 


II. 1/J EXPANSION IN THE SUPER-RADIANT 
PHASE. 

The U{l)/Z 2 Dicke model2ii2^ is described by: 

Hu{1)/Z'2 = iOaO^ a UJbJz + 

+ ( 1 ) 

where the coa , uJb are the cavity photon frequency and the 
energy difference of the two atomic levels respectively, 
the g = ^fNg^N = 2J is the collective photon-atom 
rotating wave (RW) coupling. The g' = VNg' is the 
counter-rotating wave (CRW) term. We fix their ratio 
to be 0 < g'/g = /3 < 1. If /I = 0, EqnlT] reduces to 
the 17(1) Dicke model^SiSi with the 17(1) symmetry a —>■ 


ae*®,cr“ —>■ cr“e*® leading to the conserved quantity P = 
a'I'a + Jz- The CRW g' term breaks the 17(1) to the Z2 
symmetry a —> —a, cr~ —>■ —a~ with the conserved parity 
operator H = If /; = it becomes the Z2 

Dicke modeli^ii^. 

Following^!, inside the super-radiant phase, it is con¬ 
venient to write both the photon and atom in the po¬ 
lar coordinates a = + , b = \/Xl + dpte*®'’. 

When performing the controlled 1/J expansion, we keep 
the terms to the order of ~ j, 1 and ^ 1 /j, but drop 
orders of 1/j^ or higher. We first minimize the ground 
state energy at the order j and found the saddle point 
values of Xa and A;,: 

Xa = ^^Jkl-fi^), Xb = Vj(1 - m) (2) 

where /i = uja^^b/ig + g')^- In the superradiant phase 
g + g' > gc = \/WaWb- In the normal phase g + g' < gc, 
one gets back to the trivial solution Aa = Af, = 0. 

Observe that (1) in the superradiant phase g{l -\- (3) > 
gc, A„ ~ A^ ^ j, ( 2 ) it is convenient to introduce the ± 
modes: 9± = {9a ± 9b)/2,5p± = Spa ± 6pb, A^ = A^ ± A^. 
(3) Defining the Berry phase in the -I- sector— as A^ = 
P + a where P = I, 2, • • • is the closest integer to the A^, 
so —1/2 < a < 1/2. Due to the large gap in the 0_, it 
is justified to drop the Berry phase in the — sector. (4) 
after shifting 9± —>■ 9± -|- 7 r/ 2 , we reach the Hamiltonian 
to the order of 1 /j: 

n[dp±,9±\ = ^{5p+ - a)^ + D_[Sp- +'r{Sp+ - a)f 
+ ^^aXl [sin^ 6 »_ -b sin^ 9+] (3) 

where D = ^ is the phase diffusion constant in 

the + sector, £)_ = E'jj/l&X^uJa with E'jj = (wa + Wf,)^ + 
4(g-|-g')^A^/A. The 7 = -^(1 — -^^ 7 ^^) is the coupling 
between the -I- and — sector. 

Note that the large J expansion condition is A^ ^ 1. 
In the / —^ 00 limit, it holds for any g > gc, but leads to a 
constraint on g at a finite j. In the superradiant regime, 
one can simply set sin^ 9± ^ 9^ in EqnjS] which becomes 
a quadratic theory. It can be easily diagonalized and lead 
to one low energy gapped pseudo-Goldstone mode and a 
high energy gapped optical mode. Setting /3 = I recovers 
the results for the Z2 Dicke model in the superradiant 
phase at A = 00 ini^. 

If one neglects the quantum fluctuations of the 0_ 
mode, namely, by setting 9- at its classical value 9- = 0, 
EqnOis simplified to: 

'H+ [Sp+ ,9+] = ^{Sp+-a)^ + 2ujaXl sin^ 9+ (4) 

where = ih. In the superradiant limit 

4A^Y^ ^ identify the approximate atomic 
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FIG. 3. The bound states and the quantum tunneling pro¬ 
cesses in the QT regime in Fig[2] The atomic energy ta and 
the tunneling energy Aq are shown. But the higher optical 
energy to is not shown, (a) At g = gd > gc = ^yoJaOJb/{l +P), 
the double well potential in the 9+ sector in Eqn|T] holds just 
one bound state with energy Ei denoted by a red dashed line. 
The blue dashed line shows the quantum tunneling between 
the two bound states which leads to the splitting Aq listed in 
Egn llhl The blue solid lines denote scattering states shown 
in EiglS^UJa (b) As the g increases further to g = g^ > gd, 
the lowest two scattering states in (a) also become the sec¬ 
ond bound state with energy E 2 = Ei -\- ta, the third bound 
state shows up at gr = g^s > gc 2 > gd and so on ( see also 
FigSl). The excitation energy is the ’’atomic” energy Sa- As 
explained below Egn llSI the ground state could also be odd 
parity depending on (—1)^. 


mode: 


a;% = AuJaXl 


2/3 
1 + ^ 


D = 




Ejjl + ! 


[(5 + 5 'r-Scl (5) 


which is nothing but the pseudo-Goldstone mode due to 
the CRW g' term. 

Note that for small /3 < 1, the condition to reach the 
superradiant regime ^ 1 is more stringent than 

the large J expansion condition A^ 1. 

By neglecting the quantum fluctuations of , the high 
energy optical mode in the 9- sector can not be seen any¬ 
more in Eqnl31 The approximation may not give very 
precise numbers to physical quantities, but do lead to 
correct qualitative physical picture, especially the topo¬ 
logical effects due to the Berry phase in all the physical 
quantities. 


III. 17(1) REGIME AND THE FORMATION OF 
CONSECUTIVE BOUND STATES IN QUANTUM 
TUNNELING REGIME. 

As the potential in the 9+ sector in Eqnl4] gets deeper 
and deeper, there are consecutive bound states forma¬ 
tions at Qc < Qci < gc 2 ■ ■ ■ leading to the ” atomic” en¬ 
ergy scale Cq. The QT regime in Eigl5]is signatured by the 
hrst appearance of the bound state after which there are 
consecutive appearances of more bound states at higher 
energies ( Figl5^,b ). The regime gc < g < ffci is the 
17(1) regime in Figl2] 

One can calculate all these gc < gd < 5 c 2 • • • by us¬ 
ing the Bohr-Sommerfeld quantization condition for a 
smooth potential Jl^pd9 = (n + lj2)Trh, n = 0 , 1 , 2 , • • • 


where p = and the i3„+i = 

Ei,E 2 , . is the {n + 1) — th bound state energy in 

FiglSKjb. From EqnlH we can see the m = -^,V{9) = 
“ cos 20 +) and the a and b are the two end 
points shown in FiglHK. We find that the bound states 
emerge at 


^ = (n-b l/2)^?i, n = 0,1,2,--- (6) 

In EqnlH setting n = 0, one can see when < 
jh, there is no bound state. This is the 17(1) regime in 
FiglU Substituting the expression for the phase diffusion 
constant D and the atomic mode uj-q in Eqn[S] leads to 
the condition for the 17(1) regime: 


4A 


2 

a 




(7) 


Note that for large J expansion to apply, one only need 
to require A^ > 1, So for sufficiently small g'/g = /3, 
there is an appreciable 17(1) regime gc < 9 < Pci before 
the quantum tunneling (QT) regime in Figl2] 

In this 17(1) regime, the second term in EqnH] breaks 
the 17(1) symmetry to Z 2 symmetry, the Goldstone 
mode at V = 00 simply becomes a pseudo-Goldstone 
mode^Sri^. But its effects at a finite N is much more del¬ 
icate to analyze. One can treat the second term in Eqnl4] 
perturbatively either by a non-degenerate at a 7 ^ 0 or 
degenerate perturbation expansion at a = 0. The total 
excitation P is not conserved anymore and is replaced 
by the conserved parity B = (—1)^, the energy levels 
are only grouped into even and odd parities in FigHl As 
shown in2ij at a given sector P, m = 0, ±1, • • • , ±P at 
a = 0 , a first order degenerate perturbation at m = ±1 
leads to the maximum splitting at a = 0 in Figldja: 

^m=±l.(7(l)(a - 0) - y- ^ 1+ (3 

Using Eqnl^ one can rewrite A = -j[l — |(^^^)^] > 0. 
For general ±m, one needs m — th order ( with the con¬ 
straint \m\ < P in a given sector P ) degenerate per¬ 
turbation calculation to find the the maximum splitting 
at a = 0 between the m and m + \ crossing in 

FigHt*: 


— 0) — D{m + 2 ^ ~ Rm+l — Rm (9) 

where Rm ^ = 0,1, ....P; i?o = 0 is the 

gap opening at the to -|-1 crossing in the U ( 1 ) limit Figl4^. 
Setting TO = 0 recovers EqnlHl 

Note that the degenerate pair (to, —to — 1) at the edge 
at a = ± 1/2 has different parity, so will not be mixed 
in any order of perturbations^i. It is easy to compute 
the edge gap at a = ± 1/2 by a non-degenerate pertur¬ 
bation. Obviously, the second term in EqnU] connects 
only Ato = ± 2 , so the first order perturbation vanishes. 
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one need to get to at least second order non-degenerate 
perturbation^!. There could also be a slight shift in 
the crossing point between the two opposite parities at 

(m, —771 — 1) 

Am.c/(i)(a = ±1/2) = 1) - S'(m) (10) 

where S'(m) ~ ^ leads to D,2D,3D . at 

m = 0,1, 2.at the U (1) limit /3 = 0 shown in FigH^. It 

is easy to see that at a given m — th doublet the maximum 
gap at Of = 0 is smaller than the edge gap at a = ± 1 / 2 , 
= 0) < = ±1/2). But both are of 

the same order. 

Nonetheless, the important phenomena of Goldstone 
and Higgs mode at the U{1) limit can still be observed 
in this 17(1) regime after considering these degenerate 
and non-degenerate perturbations. Various photon cor¬ 
relation functions in this 17(1) regime can be evaluated 
straightforwardly!!. 

Now we follow the formation of the bound states 
just after the 17(1) regime. When = jh, namely 
g = gel, it just holds the first bound state with Oh = 
-Oa = ^, El = 2uJaXl^. When < ^ < 
it holds the first bound state ( Figl3^ ) with energy 
El = 2 a;„Ag^sin^ 9, where ^F{0a) = ^h,F{0a) = 

d0\/ sin^ 0a — sin^ 0, 0 < 0a < 7 r/ 2 . It is easy to see 
O<F{0a)<l,F{0a = 7r/2) = l. 

I 


IV. QUANTUM TUNNELING BETWEEN THE 

TWO BOUND STATES: BERRY PHASE AND 
INSTANTONS. 

The instanton solution for a Sine-Gordon model was 
well known!!. From EqnlU we can find the classical in¬ 
stanton solution connecting the two minima from 0 + = 0 
or 0 _|_ = tt: 0+{t) = where tq is the 

center of the instanton. Its asymptotic form as r —>■ oo 
is 0+(t —>■ oo) —>• TT — The corresponding 

classical instanton action is: 

= ^ (13) 

The instanton problems in the three well known sys¬ 
tems (1) a double well potential (DWP) in a (f>'^ theory 
(2) periodic potential problem (PPP) (3) a particle on 
a circle (POC) are well documented in!!. The tunnel¬ 
ing problem in the present problem is related, but differ¬ 


When = ^h, namely g = gc 2 , it just holds the 
second bound state with 0b = —0a = f, E 2 = 2u!aX'^Y^. 
While the first bound state energy is given by F{0a) = 
1/3, El = 2uja\l-^ sin^ 0a- When ^ < ^h, it 

holds two bound states ( FiglSb ) with the energies: 

El = 2uJaXl^^sm'^ 01a, ‘^F{0ia) = 

E 2 = 2u;aXl^^Sm‘^02a, ‘^F{02a) = (H) 

where one can identify the first atomic energy in FigISb: 

28 

ea = E 2 - El = 2waA^ Y^p^(sin^ 02a - sin^ 0ia) (12) 

As expected ea is different than the uj-o in EqnlSl 
As 5 —>■ 00 , D —^ 0, while W-o —>■ so the 

left hand side of Eqnl 6 ]diverges, there are infinite number 
of bound states shown in EigH^. 

Note that because the bound state is either localized 
around 0+ = Q 01 0+ — tt, so the Berry phase a in Eqnl4] 
plays no roles, so can be dropped. However, as to be 
shown in the following section, it does play very impor¬ 
tant roles in the quantum tunneling process between the 
two bound states shown in EigI3]and[S] 


ent than all the three systems in the following important 
ways:. (1) the potential V{0) = 2waA^^^(I — cos20+) 
in Eqnl4]is a periodic potential in 0+. In this regard, it is 
different than the (j)^ theory, but similar to PPP. (2) The 
0+ is a compact angle confined in 0 < 0+ < 27r. In this 
regard, it is different than the PPP, but similar to POG. 
(3) There are two minima inside the range 0 < 0+ < 27r 
instead of just one. In this regard, it is different than 
the POG, but similar to the theory. So the present 
quantum tunneling problem is a new class one. Further¬ 
more, it is also very important to consider the effects 
of Berry phase which change the action of instanton to 
Sint = Sn + iaiT, that of anti-instanton to Sint = So — iair 
( FiglS] ) where —1/2 < a < 1/2 is the Berry phase in 
EqnSl 

Taking into account the main differences of the present 
QT problem from the DWP, PPP and POG studied 
perviously!!, especially the crucial effects of the Berry 
phase, we can evaluate the transition amplitude from 
0 + = 0 to 0+ = TT in FigO 
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FIG. 4. The atomic energy level evolution from the t/(l) regime to the QT regime as g changes at a fixed /3. The ground state 
energy has been subtracted. D is the phase diffusion constant. In the U{1) regime in (a) and (b), there are P doublets along the 
vertical ( blue dashed ) lines, (a) The energy spectrum at the U{1) limit /3 = O^i. (b) The energy level repulsions between the 
same parity at a = 0 and level crossings at a = ±1/2 between the opposite parities in the U{1) regime. The maximum splitting 
at a = 0 and the gap at the edge a = ±1/2 are given by Eqn[9] and [10] respectively. All the states are scattering states. If the 
total excitation P is even, along the vertical ( blue dashed ) lines, the doublets are organized as (e, o), (o, e), (e, o), (o, e), ■ • ■ 

( see also Fig|6^ ). If P is odd, one need flip all the parities. See Fig|7^|8^ at = 2,/3 = 0.1. If we follow the states with 
P, P + 1, P + 2.... at 1 = 0,1, 2.. respectively, then the doublets are organized as (e, o), (e, o), ■ • ■, there are relatively large shifts 
of zeros to the right delineated by the black dashed line. Only the splittings at m = ±1, ±2 at P = 1,2 are shown here. See 
FigOrlSIa at A = 2, /3 = 0.5. (c) The doublet at / = 0 becomes the first bound state ( Schrodinger cat ) denoted as B at 
g = gd, while the states at I — 1... remain scattering states as shown in Fig|2^|S^. (d) The doublet at 1 = 1 becomes the second 
bound state at y = g^, the two bound states at 1 = 0, 1 are connected by nearly straight boundaries at q = —1/2, 0,1/2. The 
states at Z = 2,.... remain scattering states as shown in FiglSja. See FigjTtlHJ; at A = 2, /I = 0.9. (e) Finally, more states become 
bound states a.t g — gcz, ■ ■ ■ where the energy level pattern becomes (e, o), (e, o), • • • in the QT regime ( see also FigjSJr ). If 
P is odd, one need also flip all the parities. There is a shift by exact one period from (b) to (e). Intuitively, a leaning tower 
with infinite stories gradually becomes straight as the system evolves from the U(l) regime to the QT regime. The maximum 
splitting at Q = 0 are given by Ean l8l9l for scattering states, Egn llOllTl for bound states ( Schrodinger cats ) and decrease from 
(b) to (e). The edge gap at a = ±1/2 are given by Egn llOl and the atomic energy in FiglJJs for scattering states and bound 
states respectively and increases from (b) to (e), but starts to become flat since the second bound state forms at Z = 1 in (d). 
The maximum gap at o = 0 is comparable to the edge gap at a = ±1/2 in the Z7(l) regime (b), but much smaller in the QT 
regime where A; <C uJa (d). This evolution from (b) to (e) is precisely observed in the ED in Fig|7]and fine structures in all 
the doublets at Z = 0,1, 2 in FigjS] The higher energy optical modes are not shown. As said in the caption of Fig|3 the U(l) 
regime disappears if /3 gets too close to 1 as shown in Fig|71:|SJ;. 


ni,n2 


ni\ 


7121 


ni —n2 ,odd 


, a;_or/2lr 2Jifre-®0 


COS aTT _ ^—2JKTe cos otti 


(14) 


where ni ( n 2 ) is sum over the instanton ( anti-instanton 
), J = (5'o/27r/i)^/^ is given by the instanton action in 
Eqnliniand K is the ratio of two relevant determinants to 
remove the zero mode of the instantons due to its center 
To listed above EqnjTSl It can be shown that K = Cquj-o 
where Co = 2 is extracted from the asymptotic form of 
the instanton solution in the r —>■ oo limit listed above 
Ean ll3l Similarly, in finding the transition amplitude 


from 0 back to 0, one only need to change 6 ni-n 2 ,odd to 
Sni-n 2 ,even in the first line, consequently, the — sign to 
the ± sign in the second line in EanlTil 


The two transition amplitudes lead to the Schrodinger 
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e^=7t 



e^=7u 



FIG. 5. (a) Scattering states: There is a global phase wan¬ 

dering around the 6+ circle from 0 to 2?! with the phase dif¬ 
fusion constant D. See Fig[4}3,Fig[8^,b. (b) Bound states ( 
Schrodinger cats ): The quantum tunneling process due to 
the instantons between the two bound states at 0-|- = 0 and 
6+ = TT. The counterclockwise ( blue dashed line ) tunneling 
is an instanton. The clockwise ( black dashed line ) tunneling 
is an anti-instanton. It is the Berry phase which leads to the 
oscillations of parities in the ground state and excited states. 
See Fig|4l:,d,e,Figl8};. 


” Cat ” state with even/odd parity: 

|e)o.sc = = 0) + l^+ = 7^)): 

|o)o,sc = -^(1^+ = 0) - l^+ = ^)): (15) 

where the overlapping coefficient is = |(a; = 0|n = 
0)1^ = They have the energy Se/o = T 

Ao/2 with the splitting between them given by: 

Ao(q:) = 8a;_o(cosQ;7r)((^(^)^/^e ^ 

ttD 

^ (cosa7r)\/iVe“^^ (16) 

where one can see that it is the Berry phase which leads 
to the oscillation of the gap and parity in EpnlTSl It 
vanishes at the two end points a = ±1/2 and reaches 
maximin at the middle a = 0. Note that the Berry phase 
a is defined2£i2i at a given sector P. So EqnlTS] has a 
background parity 11 = (—1)^. So there is an infinite 
number of oscillating parities in Ean llSI as g increases. 

Because the Berry phase effects remain in the given 
P sector, extending the results in Ref4^, we find the 
splitting in the n-th excited bound states ( n = 0,1 in 

FigEb ): 


A„(q:) 


1 ,8a;_o 
n!^ D 


)”Ao 


( 17 ) 


where n = 0,1,2,--- with the corresponding n-the 
Schrodinger ’’Cat ” state with even/odd parity and the 
energy = (n ± ^)hw_o T A„/2: 


\e)n,sc = 

|o)„,SC = “ h)fl): 


( 18 ) 


where is the n—th bound state in the left(right) 

well in FiglHKib. Putting n = 0 and projecting it to the 
coordinate space at x = 0 recovers EqnHS] ( Note that 
the projection to x = 0 does not work for n is odd ). One 
can see the higher the bound state in the FigEl the larger 
the splitting is. The main difference than the energy level 
pattern in the 17(1) regime is that all the bound states 
have the (e, o), (e, o), • • • ( or (o, e), (o, e), • • • ) pattern 
shown in Figl4ji,e and Fig|^. This important observa¬ 
tion is completely consistent with the results achieved 
from the strong coupling expansion in^ after identify¬ 
ing n ^ 1. For example, as explained ii>^, there is an 
extra oscillating sign (—1)* in Eqn.5 in22 achieved from 
the strong coupling expansion, which is crucial to recon¬ 
cile the results achieved from the two independent ap¬ 
proaches ! 


V. PHOTON, SQUEEZING AND NUMBER 
CORRELATION FUNCTIONS 


Following the procedures for the 17(1) Dicke model at 
/3 = 0 iii^i, treating the second term in EqnH] as a small 
perturbation when /3 is small, using non-degenerate per¬ 
turbation away from a = 0 and degenerate perturbation 
near a = 0 one can evaluate Photon, squeezing and num¬ 
ber correlation functions in the 17(1) regime outlined in 
FigE^. Here we focus on calculating these correlation 
fnnctions in the QT regime outlined in FigEJj. 

The above physical pictures in the QT regime inspire 
us to decompose photon and atomic operators as: 

a = aLTz, b=-bLTz ( 19 ) 


where ol = ± 1)l = are con¬ 

fined to the left well in FigE] and the = ±1 stand 
for the Left/Right quantum wells in FigEl It is the Tz 
component which contains the important Berry phase ef¬ 
fects and the quantum tunneling process between the Left 
and Right quantum well with the tunneling Hamiltonian 
Ht = ^iTx- Note that the two states Tz = ±l here is 
the two Schrodinger Cat states of the strongly interact¬ 
ing atom-photon system instead of the two levels of an 
atom az = ±1 in the original 17(1)/Z2 Dicke Hamilto¬ 
nian Eqn[TJ The two energies A/ and should appear 
in the single photon correlation function ( which con¬ 
tain the magnitude Spa, the phase 6 a and the Ising 
correlation functions as shown in FigE] ) with the cor¬ 
responding spectral weights ^ A, ^ 1 respectively. By 
using this decomposition, we will perform the 1/J ex¬ 
pansion to evaluate all the relevant photon and atomic 
correlation functions. 

Because the Z 2 symmetry is broken in either left or 
right well in FigEl one can ignore the periodicity in 0a, 9b 






m =±2 


m=+l 


T- 

- e 

, ~3D/2 


;~D/2 

e 


312 (0_o 

n=l 



1/2 (0_„ 

n=0 



(a) 


(b) 


FIG. 6. Splittings near the Berry phase a = 0 in the (a) ?7(1) 
regime in Fig|4}3 and Fig[5li. D is the diffusion constan t^°i^^ . 
m = ±1,±2,... are the magnetic quantum numbers, (b) the 
QT regime in Fig|4ji and FiglSjr. The ground state n = 0 
and the first excited state n = 1 and their splittings Aq and 
Ai due to the quantum tunnelings of instantons subject to 
the Berry phase. The blue and red transition lines can be 
mapped out by photon and photon number correlation func¬ 
tions respectively. Shown is the background parity (—1)^ is 
even (e). When (—1)^ is odd (o), one just changes even and 
odd in the figure. 


and expand the atom and photon operators as: 

{5pa? , 


aL-\a+ tXada + ^ 


— Ab -l- iXbdb + 


6pb XbOl {Spbf 

2 Xb 2 


8 X 1 


-t- 


( 20 ) 


Using the Holstein-Primakoff (HP) representation of 
the angular momentum operator = 6^6 — J, J_|_ = 
6V2J - Wb, J_ = V2J- one can evaluate the 

atomic spin correlation functions. Here, we focus on eval¬ 
uating the photon correlation functions. 

Using Eqn|31 [SJ we can calculate the phase-phase, 
density-density and density-phase correlation functions 
in the imaginary time r: 


(0a(r)0a(O)) 

{6pa{T)6pa{0)) 

((5pa(r)6»a(0)) 




2a;_o 

UJgXl 2/3 
2w_o1 +, 

= -{5ea{T)pM) 


2a;_o 
{{Spaf) = 


^gAg 2/3 
2w_o 1 + /3 

( 21 ) 


Because the 9+ is a phase confined on 0 < 0+ < 27r, 
we can define q = e*®+. For 0+ = 0,7r, we can set q = 
Tz standing for the left/right well in Fig l3] and find the 
correlation function in the I — th state: 

(e|r.(r)T,(0)|e) = e-^'- (22) 

where A; is the splitting at I — th state in Figl3]and|ni 


From the decomposition FqnlT9l [20l [21] and Fani22l 
one can evaluate the photon correlation functions: 

2/3 




■)- 


Wg 


8w_o 1 + /3 


)]e 


— Aqt 


8^ UJa 


1 + / 
2/3 


COa 


-) + —( 


2/3 


w-o 1 + /3 


)-2]e 


— (a;_o + (Ai+Ao)/2)r 


(23) 


where the first term containing the ground state splitting 
Ao has the corresponding spectral weight ^ N, while the 
second term containing the atomic energy Cg plus the 
average of the splittings at n = 0 and n = 1 in FigjBjhas 
the spectral weight ^ 1. 

Very similarly, one can evaluate the anomalous photon 
correlation functions: 


2/3 






-Aqt 


1 2^ _ )/(z£('iAb^jlp-(‘^-o + (Ai-HAo)/2)r 

8^W_oH+;9^ OJa ^ 2/3 


(24) 


where the first term containing the ground state splitting 
Ao has the same spectral weight ^ N as its counterpart 
in Fan l231 while the second term containing the atomic 
energy eg plus the average of the splittings at n = 0 and 
n = 1 in FigjBjhas a different spectral weight ~ 1, maybe 
even different sign than its counterpart in Fan l23l 

One can also compute the photon number correlation 
functions: 


(n(T)n(O)) - (nf 


^gAg 2/3 ^_(g;_n_(Ai-Ao)/2)T 

2cu’_o 1 + /3 


(25) 

where (n) = Ag is the photon number at the ground state. 
It contains the atomic energy Cg minus the difference of 
the splittings between n = 1 and n = 0 in FiglHjand has 
a spectral weight ~ N. 

So all the parameters of the cavity systems such as the 
doublet splittings Ao(a),Ai(a) and the atomic energy 
Eg are encoded in the photon normal Fan]23l and anoma¬ 
lous Green function Fan ]24l and photon number corre¬ 
lation function Ean l25l They can be measured by pho- 
toluminescence, phase sensitive homodyne and Hanbury- 
Brown-Twiss ( HBT ) type of experiments^ respectively. 


VI. COMPARISON WITH THE RESULTS 
FROM THE EXACT DIAGONIZATION AND 
THE STRONG COUPLING EXPANSION. 

When /3 1 in EpnlT] it is not convenient to perform 

the ED in the coherent basis anymore used ini^, so we did 
the ED in the orginal ( Fock ) basis. In the Fock space, 

the complete basis is |n)|j,m),n = 0,1,2,. oo,j = 

N/2,m = —j, . ,j where the n is the number of pho¬ 

tons and the |j, m) is the Dicke states. In performing 
the ED, followingi^, one has to use a truncated basis 

n = 0,1,.rig in the photon sector where the ric 100 

is the maximum photon number in the artificially trun¬ 
cated Hilbert space. As long as the low energy levels in 
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FigEland Figl8]are well below ricUJa, then the energy lev¬ 
els should be very close to the exact results without the 
truncation ( namely, sending ric ^ oo ). However, the 
ED may not be precise anymore when g gets too close 
to the upper cutoff introduced in the ED calculation as 
shown in Fig [St. 

In Fig 13 we show the ED results for the energy levels 
for = 2 at /3 = 0.1,0.5,0.9,1. It matches precisely the 
theoretically predicted energy level evolutions shown in 
FiglD At iV = 2, when (3 < j3c ^ 0.6 ( which, in fact, only 

I 


The important relation Ean ]17l takes the same form as 
Eqn.5 in^ except the absence of the extra (—1)* factor. 
As said at the end of Sec.IV, it is this absence of extra 
(—1)^ which reconciles the results achieved from the two 
independent approaches. EqnlT6] indicates that there are 
inhnite number of zeros due to the Berry phase interfer¬ 
ence effects in the instanton tunneling process in FigjS] 
EqnlT3 indicates that the positions of the zeros are inde¬ 
pendent of n. This is indeed confirmed by the ED shown 
in FiglSt for V = 2,/3 = 0.9, Z = 0,1,2 where the po- 

I 


VII. EXPERIMENTAL DETECTION OF THE 
BERRY PHASE EFFECTS OF THE INSTANTON 
TUNNELING EVENTS 

There have been extensive efforts to realize the 
Schrodinger Cat state in trapped ions^ and supercon¬ 
ducting qubit systems^. Here, the Schrodinger Cats in 
EqnllH] with n = 0,1,2... can be prepared in the QT 
regime in Figl2^, its size can be continuously tuned from 
N ^ 3 — 9, it involves all the N number of atoms ( 
qubits ) and photons strongly coupled inside the cavity 
and could have important applications in quantum infor¬ 
mation processions. 

With N ^ 10® atoms of inside a cavity^Zr— , the 
system is essentially in the thermodynamic limit of the 
Z 2 Dicke model, so the novel physical phenomena in the 
QT regime in Figl3 at finite small N explored in this 
work are hard to observe. The experimentZZ first adia- 
batically prepared the system in one of the two bound 
states in FiglS}] by applying a small Z 2 symmetry break¬ 
ing field, then turn off and on the transverse pumping 
laser in FiglT^ and observed the coherent switch with the 
frequency uJh between the two ground states in FigOb by 
an optical heterodyne detection. As emphasized in this 
work, in order to observe the Berry phase interference 
effects, one has to move away from the Z 2 limit realized 
in the experimentsZI^— , namely, 0 < ,5 < 1. This has 
been realized in the recent experiments^ which can tune 
(3 from 0 to 1. With the recent advances of manipulating 
a few atomsSiiSZ^ the number of atoms can be reduced 
to a few to a few hundreds, then the C/(l) and the QT 


weakly depends on N ), there is always a C/(l) regime 
FigllbjC before the formations of bound states in the QT 
regime in Figl3i,b. It is the Berry phase which leads to 
the parity oscillations in both regimes. However, (3 > j3c, 
the systems get to the formations of bound states directly 
in FigI3:. It is still the Berry phase which leads to the 
parity oscillations in the QT regime after the formations 
of bound states as shown in Figl33. At /3 = 1, the Berry 
phase effects and the level crossings are pushed to infinity, 
so no parity oscillations anymore in FigI3i. 


sitions of the first N = 2 zeros only depend on I very 
weakly. So between the two zeros, at I = 0,1, 2, • • •, the 
energy levels are either in the pattern (e, o), (e, o), • • • or 
(o, e), (o, e), • • • in Figidji and FiglSb. This important re¬ 
sult completely substantiates the results achieved from 
the strong coupling expansion in^Z. The fact that the 
same fantastic phenomena are reached from two inde¬ 
pendent analytic approaches, then confirmed by ED in¬ 
dicates that the results are correct, independent of the 
1/J expansion or strong coupling expansion we made. 


regime in Figl2 Figl3^,b span a large parameter regimes. 
One can first adiabatically prepare the system in the left 
or right bound states in Fig 13] with n = 0,1,..., but still 
keep the transverse pumping laser in FiglT^, then ob¬ 
serve by the optical heterodyne detectionZZ the coherent 
oscillation probability between the two bound states: 

P{a,t) = cos An{a)t (26) 

where the A„(a) is given by EanfTTl 

In circuit QED systems, there are various experi¬ 
mental set-ups such as charge, flux, phase qubits or 
qutrits, the couplings could be capacitive or inductive 
through A, V, S or the A shaped. Especially, contin¬ 
uously changing 0 < (3 < 1 has been achieved in the 
recent experimenlj^. An shown inZi^ by tuning the po¬ 
tential scattering term X^J^a^a/j and the qubit-qubit in¬ 
teraction term uJ^/j, the critical coupling in Figl3]is 
reduced to gc = \/(wq — A^)(wf, — 2u)/(l -I- f3). We ex¬ 
pect all the interesting phenomena in the 17(1) and QT 
regime in Figl3] at a finite N = 3 — 9 qubits, especially 
the dramatic Berry phase effects in both regimes can be 
observed in near future experiments. 


VIH. CONCLUSIONS AND DISCUSSIONS 

Quantum optics differs from condensed matter physics 
at least in two important ways (1) the former mainly deal 
with finite size systems, while the latter mainly deal with 
thermodynamic limit ( or edge states if there is a bulk 
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(a) 




(b) 



(d) 


FIG. 7. ED results for the energy levels a,t N — 2, P = 0.1, 0.5, 0.9, 1 for (a)-(d) respectively. For simplicity, we only show 
LOa ~ oJb case. The parity even (e) and odd (o) are indicated. There are none, one and two level crossing(s) in the normal 
regime at 1 = 0,1,2 respectively, (a) For P = 0.1, there are considerable ranges of the Uil) regime Qc < g < Qci before the 

first bound state formation at Z = 0, then it evolves to the QT regime where the bound states start to form at Qd, gc 2 , Qcs 

( which are labeled only in (d), but could be labeled in (a)-(c) also ) at Z = 0,1,2.... as g/gc increases. There is a one to 
one correspondence to Fig[4l (b) For P — 0.5, the U{1) regime becomes much smaller, (c) For P = 0.9, the Z7(l) regime 
disappears. When expanding the doublets at Z = 0, 1, ...., as g/gc increases, there are infinite energy level crossings leading to 
the oscillations of parities at the ground states at Z = 0,1, .... manifolds shown in FigO As /I —>■ 1“, all the zeros are pushed to 
infinity, (d) The Z 2 limit P = 1. There are consecutive energy level mergings at gc < gd < gd < • • • which signify the bound 
state formations at Z = 0,1, 2.... and the energy levels become approximately flat . There are no Z7(l) regime, no level crossing 

between the even and odd parity pairs. The ground state pattern is (e, o), (e, o), ■ • ■ at Z = 0, 1, 2. Only the atomic energies 

at Z = 0,1, 2.... are labeled. As g/gc —t 00 limit, they approach to lijJa from below^. 


topological order ) (2) the former mainly study pumping- 
decay non-equilibrium systems, while the latter mainly 
equilibrium systems. In this paper, we focused on the 
first feature. The combination of both features will be 
presented elsewhere^. In studying the latter, one stress 
” More is different ” as advocated by P. W. Anderson. 
Here, to study the four standard quantum optics mod¬ 
els in the former system, we take the ” Few is tricky ” 


dual point of view^Sr— which establish the connections 
between the many body physics in condensed matter sys¬ 
tems and few body problems in quantum optical systems. 
We introduced the generic Z7(l)/^2 Dicke model^i which 
incorporates all the 4 quantum optics models as its vari¬ 
ous special limits. In this paper, we investigated the new 
phenomena in this model at a finite N from the 1/J ex¬ 
pansion which is complementary and dual to the strong 
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FIG. 8. Identify U{1) and QT regimes in the results achieved from the strong coupling expansion and the ED^. In the Log 
scale, the even/odd splitting A; for N = 2 at different j3 = 0.1, 0.5, 0.9 in the row and at Z = 0,1, 2 in the column calculated by 
strong coupling expansion ( red line ) versus that by the ED ( blue line )^. The labels e and o are the parity of the ground 
states. The first even parity is in the normal regime in (a) and (b), but becomes a bound state at gd which is before the first 
zero in (c). The next one must be the odd parity state at Z = 0,1, 2. Note also the changes of the numbers in the vertical axis 
at the Log scale telling the splitting increases as Z = 0 to Z = 1, then to Z = 2 as dictated by Ean l8l9l for the scattering states in 
(a) and (b) and by Eon l 161 and 1 171 for the bound states in (c). There are also none, one and two level crossing(s) in the normal 
regime at Z = 0,1,2 respectively. Eon llBI dictates there are infinite number of zeros. The strong coupling expansion in^ to 
N—th order only gives the first N zeros. The ED gives infinite number of zeros ( but not shown here ) after the first N = 2 
zeros which can only be achieved from higher order perturbation calculations in the strong coupling expansion, (a) The strong 
coupling results^ match well with those from ED at I — 0, but not too well at Z = 1, 2 in the first N — 2 zeros at /? = 0.1 in 
the Z7(l) regime. Even so, they seem match well the envelop of the splitting at Z = 0,1, 2 ( namely, the maximum splitting at 
0 = 0). This causes no concerns because the agreement is not expected in the Z7(l) regime too close to Qc when /3 is small. 
Compare with Fig[7K, the first N = 2 zeros at Z = 0,1, 2 are due to the scattering states in the U(l) regime in Fig[4}3. (2) The 
strong coupling results^S match very well with those from ED in the first N = 2 zeros even at j3 = 0.5. The other zeros are far 
apart from the first N = 2 zeros and out of the scope in the figure. Compare with FiglTlr, the first N = 2 zeros at Z = 0,1, 2 
are still due to the scattering states in the narrow Z7(l) regime in FigUJs. Observe that the strong coupling expansion works 
very well and reproduces precisely the Berry phase effects in the Z7(l) regime not too close to Qc- As shown in Fig|4]3, there are 
shifts of zeros to the right if one follow the ground state with the odd parity. (3) The match is essentially perfect in the first 
N = 2 zeros at /3 = 0.9 in the QT regime. The other zeros are far apart from the first N = 2 zeros and out of the scope in the 
figure. At too strong couplings, the ED may become ( noise ) un-reliable due to the cutoff introduced in the ED. Compare with 
FiglZb, all the zeros are well beyond gci,gc 2 , gcs... in Fig[7J:, so the first A = 2 zeros are the first two bound states at Z = 0,1, 2 
in the QT regime in Fig[2]and[4^. The U{1) regime is squeezed out. Their locations are nearly independent of Z as dictated by 
EanfTTl for the bound states satisfying G = S> 1. 
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coupling expansion used in^. 

It is constructive to compare the two analytic methods. 
The instanton method employed in this paper which is in 
the spirit of path integral can map out the physical pic¬ 
ture clearly. It starts from the U{1) limit with /3 = 0 and 
animate the consecutive formation of the bound states, 
quantum tunneling processes subject to the Berry phase 
effects shown in Fig ]2l3l5l the energy level evolution from 
the U{1) to the QT regime in FiglU It is the Berry phase 
interference effects which lead to the infinite oscillations 
in the parity of the ground state doublets in EqnlT6] and 
also excited state doublets EqnlT3 It can be used to 
predict the zeros happen at a = ±1/2 and the max¬ 
imum splittings happen at a = 0 phenomenologically, 
but can not used to predict where the zeros and maxi¬ 
mum splittings happen in g/gc- For example, it is hard 
to determine the behaviors of these zeros as /3 —>■ 1“ 
limit. Only when taking the results achieved from the 
strong coupling expansion from the (3 = IZ 2 limit, one 
can see that all the zeros are pushed into infinity in the Z 2 
limit. However, we need to evaluate the photon correla¬ 
tions functions Ean i23l [24l and [25] separately in the C/(l) 
regime ( FiglB^ ) by the perturbation theory and in the 
QT regime ( FiglS}] ) in an intuitive and phenomenolog¬ 
ical way. The strong coupling expansion employed in^S 
which is in the spirit of canonical quantization can not 
distinguish the differences between the scattering states 
and the bound states, therefore not the physical process 
of the bound state formation in Fig ]2l3l5l It starts from 
the Z 2 limit with 13 = 1. The Berry phase effects are 
only implicitly embedded in the expansion in term of the 
anisotropic parameter away from the Z 2 limit (3^1. So 
the physical picture is less clear. However, it can be used 
to evaluate the first N zeros very precisely when com¬ 
pared with the ED in FiglT] and FigE] It can also be 
used to calculate all the photon correlation functions in 
both the QT and C/(l) regimes systematically and in a 
unified scheme. So the two analytical methods are com¬ 
plementary and dual to each other. Their combination 
leads to rather complete understandings of both physical 
mechanisms and quantitative values of all the experimen¬ 
tal measurable quantities in the QT regimes and the C/(l) 
regime not too close to the QCP at IV = 00 in Figl2]and 

H 

At the U{1) limit /3 = 0, there are infinite level cross¬ 
ings due to the Berry phase effects at a finite N ( FigEJi 
) as presented Turning on a small /3 will only 

lead to level repulsions between the same parity states, 
the Berry phase still leads to the level crossings between 
the even and odd state, therefore the alternating par¬ 
ities on the ground state and also all the doublets at 
\m\ = 1, • • • P in the C/(l) regime ( Figldb ). When (3 gets 
bigger, the system evolves into the QT regime where the 
Berry phase continue to play a crucial role leading to in¬ 
terference between different instanton tunneling events ( 
Fig lit and d ). As /3 gets close to 1, the 17(1) regimes dis¬ 
appears, the normal state directly gets to the QT regime, 
the Berry phase effects show up after the formations of 


all the bound states. At the Z 2 limit /I = 1, A = 0, 
the Berry phase effects are pushed into infinity, so there 
is no level crossings between opposite parities anymore 
in FiglTJl, the energy levels statistics at a given parity 
sector satisfy the Wigner-Dyson distribution in the su¬ 
perradiant regime!^. However, at any /3 < 1 in Eqnlll as 
shown iit^, it is the extra term A(a^ — a)iJy which in¬ 
troduces frustrations, therefore Berry phase effects into 
the ^ 2 / 17 ( 1 ) Dicke model. They leads to infinite level 
crossings with alternating even and odd parity in the 
ground state and all the doublets at 1 > 1 ( Figldji ). 
Combining the physical picture from /3 small achieved 
from 17(1) limit by instanton method in this paper to 
large /3 ^ 1 achieved from Z 2 limit by the strong cou¬ 
pling expansion method in^Sj we conclude that it is the 
Berry phase effects which lead to the level crossings at 
any 0 < /3 < 1 except at the Z 2 limit 13 = 1. From FiglU 
we expect that the level statistics at a given parity sector 
still satisfies the Possion statistics in the normal regime, 
the Wigner-Dyson distribution in the QT regime, but 
it remains interesting to see how it changes in the 17(1) 
regime in Figl2l It was shown in Ref4i that it is Berry 
phase effects in the instanton tunneling events in the 2 ± 1 
compact QED which leads to the Valence bond order in 
2d quantum Anti-ferromagnet. Here we showed that it 
is Berry phase effects in the 0 ± 1 dimensional instanton 
tunneling events in the compact phase of photons which 
leads to the infinite level crossing with alternating parity 
in the ground and low energy excited states ( FigEE). 

There are some illuminating duality in both the Hamil¬ 
tonian and the quantum numbers characterizing the en¬ 
ergy spectrum of the 17(1)7^2 Dicke model. In the 
present paper, we start from the Hamiltonian in its 
17(1)7^2 representation EqnlUand use its complete eigen¬ 
states |l)m|w) when j3 is not too large ( namely near 
the 17(1) limit ). In the 17(1) crossover regime in FigE 
and E the Landau level index I = 0,1, ...,V ( N + 1 
Landau levels ) denotes the high energy Higgs type of 
excitation, the magnetic number m = —P,—P ± 1,.... 

( no upper bounds ) denotes the low energy pseudo- 
Goldstone mode^l. However, in the strong coupling ex¬ 
pansion used in^, we start from the Hamiltonian in its 
dual ^ 2 / 17 ( 1 ) representation and use its complete eigen¬ 
states |l)m|j, m) when 1 —/3 is not too large ( namely near 
the Z 2 limit ). In the QT regime in FigEandE the Lan¬ 
dau level index I = 0,1,... ( no upper bounds ) denotes 
the low energy atomic excitation, the magnetic number 
m = —j, —j + 1,...., j ( also 2j ± 1 = V ± 1 ) denotes 
the high energy optical mode. So the Landau level index 
and the magnetic number exchanges their roles from the 
17(1) to the QT regime. The crossover between the two 
basis is precisely described in FigE Of course, when (3 
gets too close to 1, the 17(1) regime disappears and so 
does the duality relations. 
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